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Abstract: We present that by combining Crossing Statistic [1,2] and Smoothing 
method [3-5] one can reconstruct the expansion history of the universe with a very 
high precision without considering any prior on the cosmological quantities such as 
the equation of state of dark energy. We show that the presented method performs 
very weU in reconstruction of the expansion history of the universe independent of 
the underlying models and it works well even for non-trivial dark energy models 
with fast or slow changes in the equation of state of dark energy. Accuracy of 
the reconstructed quantities along with independence of the method to any prior 
or assumption gives the proposed method advantages to the other non-parametric 
methods proposed before in the literature. Applying on the Union 2.1 supernovae 
combined with WiggleZ BAO data we present the reconstructed results and test the 
consistency of the two data sets in a model independent manner. Results show that 
latest available supernovae and BAO data are in good agreement with each other 
and spatially flat ACDM model is in concordance with the current data. 
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1. Introduction 

Reconstructing the expansion history of the universe and properties of dark energy 
is one of the main goals of today's cosmology to understand our universe and its 
components. There have been many approaches in last decade proposed to do the 
reconstruction of the expansion history and one can generalize them in two cate- 
gories of parametric and non-parametric methods. Parametric methods are viable 
approaches if we know the actual class-form of the phenomena we are studying and 
we can use them to put constraints on the parameters of the model. See [6-10] for 
details of data analysis and methods of parametric reconstruction of the properties 
of dark energy using supernovae data. However dealing with phenomena that we 
have no clear idea about its nature and behavior, using parametric methods can be 
misleading since the underlying actual model might not be covered by the assumed 
parametric form. Dealing with uncertainties in the dispersion of the data adds an- 
other complication to the analysis and leave us with no clear way to find this fact 
that we might have chosen an inappropriate parametric form. This raises the im- 
portance of the non-parametric and model independent approaches to find out the 
behavior of the phenomena in a more direct way by avoiding parametrizing cosmo- 
logical quantities [3-5,11-21]. However non-parametric approaches have their own 
shortcomings. For instance, estimation of the errors can be a tricky task in many 
cases since in some methods one cannot easily assign the degree of freedom in the 
likelihood analysis. For a review over this subject look at [22]. 

Smoothing method was proposed to reconstruct the expansion history of the uni- 
verse in a completely model independent way as a top-down approach starting with 
the data to reconstruct the cosmological quantities step by step from the luminosity 



distance dL{z) to Hubble parameter h{z) and the deceleration parameter q{z). This 
method has been used broadly since its proposition and has passed different tests 
and has shown its strength. Nevertheless, while smoothing method is a promising 
approach to reconstruct the underlying expansion history of the universe and prob- 
ably the best way to find a feature in a dataset, estimation of the error-bars on the 
reconstructed quantities have been an issue since it is not possible to define degrees 
of freedom in this approach. Error amplification because of taking the derivatives 
(moving from smooth dL{z) to h{z) and then q{z)) and at the same time controlling 
the bias in the reconstruction are another important related issues in the smoothing 
method. In this paper we use the idea of Bayesian interpretation of Crossing Statis- 
tic [1,2] to estimate the uncertainties of the reconstructed cosmological quantities 
in a well defined and robust statistical manner. Bayesian interpretation of Crossing 
Statistic [2] was introduced recently for the purpose of model selection and falsifying 
cosmological models in a purely model independent manner while we have no infor- 
mation about the underlying model of the universe. In this paper by combining the 
two methods, error-sensitive smoothing method and the Crossing Statistic, we re- 
construct the expansion history of the universe and estimate the redshift evolution of 
some cosmological quantities such as deceleration parameter q{z) and Om diagnostic 
Om{z). 

In the following we will explain first briefiy about the smoothing method and 
how it can be used for the reconstruction of the properties of dark energy. Next 
we explain about the Crossing Statistic and its Bayesian interpretation and then 
we introduce how we can combine the two methods to reconstruct the expansion 
history of the universe and set the confidence limits. Efficiency of the method will be 
tested using simulated data and it will be shown how precisely it can reconstruct the 
expansion history of the universe. In fact it is not unfair to claim that the proposed 
method works with considerably higher precision in comparison to other available 
non-parametric methods while here we are not even using any priors in the analysis. 
Then I apply the method on the current supernovae and Baryon Acoustic Oscillation 
(BAO) data. At the end I summarize and highlight the advantages of this approach 
to other available methods. 

2. Method and Analysis 
2.1 Smoothing Method 

Smoothing method is a completely model independent approach to derive the dL{z) 
relation directly from the data, without any assumptions other than the introduction 
of a smoothing scale. The only parameter used in the smoothing method is the 
smoothing width A, which is constrained only by the quality and quantity of the 
data, and has nothing to do with any cosmological model. The smoothing method 
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is an iterative procedure with each iteration typically giving a better fit to the data. 
It has been shown in [3-5] that the final reconstructed results are independent of the 
assumed initial guess, dL{ziy below. 

The modified smoothing method (error-sensitive) can be summarized by the fol- 
lowing equation [5]: 



\ndL{z,Ay = In diizY 
+N{z) - 
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where dL^z) is the data, N{z) is the normahzation factor, dL^ZiY is the initial 
guess model and A is the width of smoothing. 

The absolute brightness of the supernovae is degenerate with since the ob- 
served quantity is the distance modulus /i(^). The outcome of the smoothing method 
is therefore H^dL{z)/c = d'^£^{z) = (1 + z)D{z). In this paper I choose A = 0.30 
which is similar to the value being used in [5] . Complete explanation of the relations 
between the A, the number of data points, quality of the data and the reconstructed 
results can be found in [3,4]. It has been shown before that smoothing method is 
a promising approach to reconstruct the expansion history of the universe however, 
setting the confidence limits have been an issue and previously bootstrap approach 
was being used to set the confidence limits. In this paper, the reconstructed form of 
the dL{z) will be used as a mean function in the full reconstruction process which in- 
cludes the idea of Bayesian interpretation of Crossing Statistic as it will be explained 
in the next section. 

2.2 Crossing Statistic and Reconstructing Dark Energy 

The main idea behind this work lies on this fact that the actual model of the universe 
and the reconstructed models using smoothing method would have one or two mild 
crossings (the distance modules of the actual model of the universe and the recon- 
structed ii{z) from smoothing method would cross each other at one or two points 
in the data redshift range) and Crossing functions multiplied by the mean functions 
generated by the smoothing method would cover the actual model of the universe 
to an indistinguishable level. It is in fact combining a non-parametric method with 
a parametric method to define and set the confidence limits. Crossing function is 
defined by Chebichev polynomials [2]: 



T/,(Ci, C2, ^) = 1 + Ci(^) + C2[2{^f - 1], 

^max ^max 



(2.2) 
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and we fit l^^^^ooth^^) ^ I^Smooth{z) x T//((7i, (^2, ^) to the data and find the 
best fit point (7^^^^ C^^^'^ in the hyperparameter space and also Ci, C2 points related 
to the 1(7, 2(7 and 3(7 confident limits. Each T//(Ci,C2,^) (where Ci,C2 pairs are 
within hyperparameter confidence contours) multiplied to f^smoothi^) represents a 
reconstruction of the expansion history of the universe consistent to the data. 

We should note that though the reconstructed luminosity distances are model 
independent, the derived uncertainties do depend on the number of Chebichev terms 
included. In fact depends on the number of Chebichev terms, the Ax^ (with respect 
to the best fit point) for different confidence limits can vary and this can affect the 
width of the error-band. In this paper to estimate the errors we have included the 
Chebichev polynomials up to the second order and this selection is based on an 
argument we have discussed in [2]. The fact is that within the range of the data 
up to redshift of around z = 2, we cannot expect to have more than two crossings 
between different assumed models. This is even valid for models with significant 
differences, e.g SCDM and LCDM models that though they are so different in their 
predictions for the expansion history of the universe they do not cross each other 
in more than two points (for different values of matter density). This is due to the 
fact that the observables, fi{z) are integrated quantities over the inverse of Hubble 
parameter (expansion history) and they monotonically increases by redshift. In [2] 
we have discussed this issue in more details that the data is not sensitive to more 
than two crossing between models (equivalently we can include up to the second order 
of Chebichev polynomials for error estimation). However, this is true that in some 
cases even assuming one crossing term (using Chebichev polynomial of the first term 
only) would be sufficient for the purpose of error estimation but it is always safer to 
consider all unexpected possibilities. In other words, having a broader error-bands 
which contains the fiducial model with a very high certainty would be more reliable 
and acceptable than having a very narrow error-band for the reconstructed results 
but with possibilities of missing the fiducial model at some redshift ranges. As we 
will see in the results, in all assumed cases for different models of dark energy (even 
for the non trivial kink model) the reconstruction method works pretty well and the 
fiducial models are covered by the error-bands at all redshift ranges. 

To test the method we apply first the procedure to the simulated data. I assume 
three sets of simulated data based on three dark energy models (all spatially fiat 
universes). First model is an evolving dark energy model or Kink model with the 
same parameters used in [2,24]. Equation of state of dark energy in this particular 
model is given by: 



f\ ^( , 1 + exp A,^ 1 + 2;, ^ 

W[Z) = ^0 + (^m - ^0) (2.3) 

1 - exp(A^ ') 
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X 1 



exp(A,-')+exp(A,-'(l + ^t)"') 



with the constants having the values Wq —1.0, —0.5, Zt 0.5, 0.05. 

I assume = 0.27 in the simulations. Second model is based on a dark energy 
model with constant equation of state of w{z) = 0.9 and with = 0.27 and finally 
the third simulated data is based on ACDM model with w{z) = —1 and = 0.27. 
Simulated data are based on future space based supernovae data with 2298 data 
points in the range of 0.015 < z < 1.7 with intrinsic dispersion of = 0.13 [25]. 
The Kink model I used in this analysis has a special form of the equation of state of 
dark energy that at low redshifts it converges to w{z) = —1 and at higher redshifts 
it smoothly changes to w = —0.5. This evolving equation of state of dark energy is 
not reducible to the common form of w{z) = Wq + Waz/{1 + z) [26] or many other 
dark energy parametrizations. Hence using usual parameterizations results to wrong 
reconstruction of dark energy equation of state if the data is based on the Kink model 
(look at upper-right panel of figure. 6 in [24]). 

Deriving dL{z) we can accordingly derive h{z)^ Om{z) and q{z): 



where derivatives are respect to redshift and h{z) = H{z)/ H^. To derive all these 
quantities we do not need to know the value of the matter density and in particular 
Om{z) can be used directly to falsify the cosmological constant which requires Om{z) 
should be a constant at all redshifts [15]. 

We should note that the observables that are tightly related to the lumi- 
nosity distances dL{z)^ are integrated properties that are functions of the evolution 
of the universe and there might be diflFerent dark energy models (assuming different 
matter density and curvature) that result to the same observables. The concept of 
cosmographic degeneracy does exist and make it almost impossible to distinguish 
between some certain dark energy models if we have no tight constraints on the cur- 
vature and matter density [27]. However, in this paper we try to avoid this issue 
and we work on h{z) and q{z) directly. The reason is that unlike the equation of 
state of dark energy w{z)^ these quantities h{z) and q{z) are direct derivatives of 
the reconstructed distance results and they are independent of the actual value of 
the matter density. We have also assumed a fiat spatial curvature through out the 
paper. 

To have an idea of the reconstruction process, in figure. 1 we show a set of 
simulated data and the reconstructed /x(z) and its error bands. The simulated data 
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Figure 1: A set of simulated data and the reconstructed and its error bands. The 
simulated data is based on a spatially flat LCDM model and the actual model is subtracted 
from the data and the reconstructed results for better clarification. 



is based on a spatially flat LCDM model and the fiducial model is subtracted from 
the data and the reconstructed results for better clarification. In this paper by la, 
2a and 3cr we mean 68%, 95% and 99% confidence limits respectively. 

In figure. 2 we see the reconstructed h{z)^ Om{z) and q{z) for the three assumed 
models. One can clearly see that the reconstruction procedure performs perfectly 
well in reconstructing the cosmological quantities. One should notice that here we 
are not assuming any particular prior on any cosmological quantity which gives it an 
advantage to other reconstruction processes. 

2.3 Real Current SN la and BAO Data 

Now we can apply the method on the current available data. We use Union2.1 
supernovae compilation [28] that is the most recent compilation of the supernovae 
data along with WiggleZ baryon acoustic oscillation data [29] . Using BAO data we 
consider the dz ratios all scaled to rf^ at z = 0.106 to make the BAO data independent 
of the knowledge of the early universe (look at table. 1). dz is given by: 

4 - (2.7) 

where rs{zcMB) is the sound horizon at the epoch when CMB photons decouple 
from baryons and Dy{z) is the dilation-scale distance. 



Dv{z) 



dhi^Y ( cz 



1/3 
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_(l + z)2 \H{z)^ 

In fact by using the dz ratios we remove the rs{zcMB) from the equations and 
minimize the effects of the assumptions of the early universe from our analysis. 
Similar to the procedure we used in the previous section we apply first the smoothing 
method on the supernovae data to reconstruct the mean functions and then we use 
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Figure 2: Reconstructed h{z), Om{z) and q{z) with their la (green), 2cr (blue) and 3cr 
(magenta) confidence hmits for three sets of the simulated JDEM data based on a kink 
dark energy model on the top, ACDM in the middle and w — —0.9 in the bottom, all 
with fiom — 0.27. Strength of the method in reconstructing the cosmological quantities is 
obvious. Black double-dot lines represents spatially flat ACDM model with fiom = 0.27 for 
comparison. 



these mean functions along with Crossing functions fitting supernovae and BAO 
data to find out the Crossing hyper parameters and define the confidence limits. 
Reconstructed Om{z) and q{z) using Union 2.1 data is shown in figure. 3 upper 
panels. In the lower panels we can see the same quantities using combination of 
Union 2.1 supernovae and WiggleZ BAO data. One can clearly see that the standard 
flat ACDM model has a proper consistency to the data. In figure. 4 we see the la 
and 2a confidence contours of Ci , C2 hyperparameters using supernovae data and 
BAO data individually. One should remember that these are the coefficients of the 
Crossing function around the mean function generated by smoothing method using 
supernovae data. So it is expected that the Ci , C2 contours using supernovae data be 
centered around 0,0 point. This plot shows that the BAO and supernovae data are 
in proper consistency with each other without assuming any cosmological model. It 
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Inverse Covariance Matrix 









dz=0.10Q 
dz=0.20 


dz=0.106 
dz=0.35 


dz=0.106 
dz=0.44 


dz=0.106 
dz^O.60 


dz=0.106 
dz = 0.73 


Mean Values 




=0.106/ dz^ 


=0.20 


284.709 


-103.845 


-4.773 


-9.600 


-6.961 


1.764 ±0.097 


dz-- 


=0.106/ dz^ 


=0.35 


-103.845 


91.940 


-2.552 


-5.133 


-3.722 


3.063 ±0.170 


dz-- 


=0.106/ dz^ 


=0.44 


-4.773 


-2.552 


14.546 


-9.575 


2.275 


3.668 ± 0.328 


dz-- 


=0.106/ dz^ 


=0.60 


-9.600 


-5.133 


-9.575 


30.340 


-10.771 


4.628 ± 0.299 


dz-- 


=0.106/ dz^ 


=0.73 


-6.961 


-3.722 


2.275 


-10.771 


12.954 


5.676 ± 0.398 



Table 1: WiggleZ baryon acoustic oscillation distance ratios scaled to dz=o.i06 ^^nd their 
inverse covariance matrix [30]. Using the BAO distance ratios one can reduce the effects 
from the early universe on the distance measurements. 



also shows that supernovae data is still a much more powerful probe of the expansion 
history of the universe in comparison to BAO data (notice the size of the confidence 
contours). One reason lies on the fact that supernovae data covers a very large 
redshift range (from as low as z = 0.015 up to about z = 1.5) and also the number 
of data points are considerably larger in the supernovae compilations. 

3. Conclusion 

In this paper we present a robust and easy to use method of model independent 
reconstruction of the expansion history of the universe which performs perfectly well 
for any dark energy model. Combining the smoothing method and Crossings Statis- 
tic results in an approach which can be easily used to reconstruct the expansion 
history of the universe and the properties of dark energy without setting any prior 
on the cosmological quantities. This is an important advantage over other avail- 
able reconstruction methods in which one has to set some priors on the equation of 
state of dark energy or some other cosmological parameters. The presented approach 
also resolves the problem of defining the confidence limits around the reconstructed 
cosmological quantities using smoothing method. It has been shown before that 
smoothing method is a very powerful approach to reconstruct the expansion history 
of the universe but because of its non-parametric nature, it was not possible to de- 
fine the confidence limits in a straightforward way. In this work, using the idea of 
Bayesian interpretation of Crossing Statistic this problem is resolved and one can 
set proper confidence limits on the reconstructed quantities. We have shown the 
strength of the method in reconstruction of the expansion history of the universe 
using simulated data where we could reconstruct different dark energy models, even 
some non trivial ones, with a very high precision. We have also used the method to 
test the consistency of the two important cosmological data sets. Union 2.1 super- 
novae data and WiggleZ BAO data in a completely model independent way. These 
two data sets are pretty well consistent with each other, however, at the current 
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Figure 3: Reconstructed Om{z) and q{z) with their la (green), 2(j (blue) and 3cr (ma- 
genta) confidence limits using Union2.1 supernovae data (on the top panels) and using 
combination of Union2.1 supernovae and WiggleZ BAO data (on the bottom panels). One 
can see that including BAO data into the analysis does not change the reconstructed re- 
sults significantly. At the current status of the cosmological observations supernovae data 
contain much more information on the expansion history of the universe in comparison to 
BAO data. 

status of cosmological observations supernovae data has much more power in recon- 
struction of the expansion history of the universe. Based on our analysis spatially 
flat ACDM model has a proper concordance to the given Union 2.1 supernovae and 
WiggleZ BAO data, even though, we cannot still rule out many other dark energy 
models including other two models discussed in this paper. 
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